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TITLE: 

TYPES 
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TEMPORARY STORAGE; 
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DURATIONS 

ACCURACY; 
DESCRIPTION; 



NOTE; 



METHODS 



Regular Spherical Bessel Functions 

Closed 

159 (12^ plus S3 and .T5O 

0-7 plus N+l locations at S3 

S3 contains the address where the table of functions 

j Q (x), j-^x), ,.., j n (x) is to be placed. 

Variable (approx * 200 ms . ) depends on N and the relation 

ship between x and N 

Maximum error 2 

Enter with scaled argument 2 x in the accumulator and 

the program parameter; 



10 « Usually between 10 " and 10 



p + 1 



50 N F 
50 PF 



26 qF 



fe < X < 19 
N < 31 



where q is the location of this program and N is the 

maximum order of the Bessel Function to be found, N 

should not be less than 3. If x is greater than N + l/2, 

N will be changed to the next integer larger than x . On 

exit j-u.(x) will be found in kS3» 

This program contains Logrithm routine S3 at 12^L and 

Sine-Cosine routine T5 at I38L* These are available 

for use in other connections. Merely enter at 12^L for 

S3 and 138L for T5« 

A complete discussion of the method may be found in the 

"Journal of the Association for Computing Machinery" 

Vol. 6, page 366, (July, 1959)* 

This ppogram starts by choosing a k greater than N and 

assumes that J = and J, , is arbitrary. It then 

computes the ratios — by using the 



r = 
n 



'n+l 



n 



- 2 - 

recursion relation r n = — — — ' ' 5 r n = 0« 

2n + 1 » xr 
n 

When r becomes greater than one,, the program starts 

to use the recursion relation for j's; 

t 2n + 1 - t 

J n-1 " x 3 n " J n+1 

The j's so found will be proportional to the true j s s* 
the proportionality constant <X is 

oC - -— = (j Q ~ x^) cos x + xj Q sin x. 

J n 

The true j's are then found using oL and the ratios 
found before • 

The ratio r become^ greater than one when the function 

n ° 

leaves the monotonic region and begins to oscillate . 

Therefore r is not used in the oscillating region since 
n 00 

r could blow up at the zeros of the function* 

n 

One should note that r = r « 

n n 

ERRORS % This method is not subject to the large errors to which 

the simpler method of upward recursion from j n and j, 
is subject* Since the ratios are always of the order 
of unity., one retains full accuracy* The j's and j's 
can be made as proportional as one wishes by choosing 
k to be sufficiently high* 
This program was checked using existing tables of the 



functions* For x less than 10, the errors were about 
10 , For x greater the 
accuracy than the tables , 



10 • For x greater than 10, this program has greater 



DAT E June 9, 1961 
PROGRAMMED B Y R„ Parsons 

APPROVED md^A/J^yU^^j^ 



HJ 



ON THE FOLLOWING FIVE- PAGES ARE TABULATED REPRESENTATIVE VALUES FOR THE 
REGULAR SPHERICAL BESSEL FUNCTIONS, THE LEFT COLUMN GIVES VALUES COMPUTED 
BY ILLIACc THE RIGHT COLUMN ARE THE VALUES GIVEN IN TABLES OF SPHERICAL 
BESSEL FUNCTIONS (VOLUME l). IN NO CASE IS THE ERROR GREATER THAN 2-10" 1:L o 
IN ALL CASES WHERE THE TABLES DO NOT GIVE THE VALUE TO ELEVEN FIGURES 9 THE 
ILLIAC VALUE ROUNDED IS THE VALUE GIVEN IN THE TABLE. 

THE VALUES ARE TABULATED FOR A GIVEN x IN THE FOLLOWING ORDERS J (x), ^(x 
*„. 5 j 15 (x). VOLUME 1 OF THE TABLES DOES NOT LIST ORDERS HIGHER THAN 15 . 



4 - 



X = ol 



ILLIAC 

+998334166^76 

+033300011894 

+000666190606 

+000009518519 

+000000105771 

+OOOOOOOOO96O 

+000000000007 

+000000000000 

+000000000000 

+000000000000 

+000000000000 

+000000000000. 

+000000000000: 

+000000000000 

+000000000000 

+000000000000 



TABLES 

+99833^1665 

+03330001190 

+000666190606 

+0000095185197 
+000000105772 
+000000000961 
+00000000007^ 

+000000000000 
+000Q00000000 
+000000000000 
+000000000000 
+000000000000 
+000000000000 
+000000000000 



x = .2 



ILLIAC 

+993346654015 
+066400380676 

+002659056079 
+000076021317 
+000001696045 

+000000030735 

+000000000473 

+000000000005 

+000000000000 
+000000000000 
+000000000000 
+000000000000 
+000000000000 
+000000000000 
+000000000000 
+000000000000 



TABLES 

4.9933466540 

+06640038067 

+002659056080 

+000076021318 

+000001690046 

+000000030756 

+000000000473 

+000000000006 

+000000000000 
+000000000000 
+000000000000 
+000000000000 
+Q00000000000 
+000000000000 
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x = «5 



ILLIAC 

+958851077232 
+16253705^3^ 
+OI637HO6605 
+001174035^2 
+OOOO65389606 
+000002977^66 
+000000114664 

+000000003825 
+000000000111 
+000000000002 
+000000000000 
+000000000000 
+000000000000 
+000000000000 
+000000000000 
+000000000000 



TABLES 

+9588510772 

+1625370306 

+01637110661 

+001174035444 

+000065389606 

+000002977^67 
+000000114665 
+000000003826 

+000000000112 

+000000000003 

+000000000000 
+000000000000 
+000000000000 
+000000000000 



x = 1,0 



ILLIAC 

+84i470984886 
+301168678945 
+0620350520Q9 
+0090065811^6 
+001011015807 
+000092561158 
+000007156936 

+000000^79013 

+000000028263 
+000000001490 
+000000000071 

+000000000002 
+000000000000 
+000000000000 
+000000000000 
+000000000000 



TABLES 

+84l4709848 

+3011686789 

+06203505201 

+00900658112 

+001011015808 

+000092561158 

+000007156936 

+000000479013 

+000000028265 

+000000001491 

+000000000071 

+000000000003 

+000000000000 
+000000000000 



x = 3.0 



ILLIAC 

+0^704000272o 
+3^5677^997^8 

+298637^97071 
+152051662026 
+0^614971^326 

+016397^80955 
+00397^382509 

+000824843251 
+000149833755 
+000024214698 
+000003526004 

+000000467327 
+OOOOOOO56845 
+OOOOOOOO639O 
+000000000668 
+000000000064 



TABLES 

+0470400027 
+345677^998 

+298637^971 
+1520516620 

+0561497L433 

+01639748096 

+00397^382510 

+000^484^3253 

+000149833756 

+000024214699 

+000003526004 

+000000467328 

+000000056846 

+000000006390 



x = 5*0 



ILLIAC 

-191784854944 
-095089408076 

+134731210115; 
+229820618220 
+187017655388 
+1068lll6l479 
+047966899867 
+017902778180 
+005741434674 

+00161809971^ 
+000407344243 

+0000927^6110 
+OOOOI9287863 
+000003693207 
+OOOOOO655453 
+000000108426 



TABLES 

-1917848549 
-0950894081 
+1347312101 
+2298206182 
+1870176553 

+io68ni6i46 
+0^796689986 
+01790277818 
+005741434675 

+001618099715 
+000407344244 
+0000927^6110 
+000019287863 
+000003693207 
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x = 10.0 



ILLIAC 

-05^02111078 
+078^66941812 
.^0779^2193615 
-039^9584^992 
-105589285105 
-05553^511603 
+0^5013223^9 
+1133862 3065f 

+125578023653 
+100096^09551 
+06^0515^49^ 

+017215999744 
+007465584475 
+00294l078|4l 
+001063542713 



TABLES 

-05440211109 

+07846694 

+07794219363 

-03949584 

-1055892851 

-05553451162 

+04450132233 

+11338623 

+1255780356 

+10009646955 

+06460515449 

+03557441 

+017216000 

+007465584 



x = 15.0 



ILLIAC 

+043352522667 
+053536029051 
-032645316858 
-064417801330 
+0025836762^4 
+065968007077 
+045792862291 
-0262808591% 
=0^2073722049 

-O55402691882 
+OOI896979OO5 
+058058462489 
+O87125996813 
+087151532201 
+O6974676115O 
+047692206019 



TABLES 

+04335252 
+05353603 
-03264532 
-0644l780 
+OO258368 
+O65968OI 
+04579286 
-O2628086 
-O72O7572 
-05540269 
+OOI89698 
+05805846 
+08712600 
+08715153 
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x = 19.O 




ILLIAC 




TABLES 


+00788827^196 




+00788827 


-0516219128^2 




-05162191 


-016039102537 




-01603910 


+0^01096386 




+0^7^110 


+03350266^361 




+03350266 


-031531^13262 




-031531^1 


-O51757693096 




-05175769 


-00388l7 i i-5173 




-00388175 


+048693157^36 




+04869316 


+0^^9307092 




40^71^931 


-0012i|-38503 I i-5 




-00124385 


-Olf8824o89050 




-0^882409 


«D578589^te9l 




-O5785899 


-0273P6166596 




-02730617 


+019055^9^336 






+056390868553 







LOCATION 



ORDER 



NOTES 
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3 

k 

5 
6 

7 
8 

9 

10 
11 
12 

13 
Ik 

15 
16 



OOK(Vl^) 
40 3F 
iH kF 
K5 F 
^2 110L 
10 20F 
k2 kF 
50 111L 
L5 ^F 
00 IF 
kO F 
L5 3F 
10 33F 
FO F 
32 7L 
22 9L 
L5 3F 
10 3^F 
Fk 111L 
kO kF 
50 111L 
L5 ^F 
00 IF 
Fk 111L 
00 33F 
kO F 
L5 3F 

66 F 

S5 S3 
^0 F 
50 F 
7J F 
kO IF 
10 IF 
Lk 112L 




2" 5 x — -> 3 



Plant link 



N .• 2~ 39 — ) k 
q = 



N(0) = n . 2' 



■38 



x - (n + 1/2) • 2 
> 
< 



n = £x] + 1 



-38 



(N + 1/2) • 2~ 5 



N(0) = 



1 x 

U s= 



n + 1/2 



N(l) = u 



LOCATION 



ORDER 



17 
18 

19 
20 
21 
22 

23 
2k 

25 
26 
27 
28 
29 
30 

31 
32 
33 
3^ 



kO 2F 
50 2F 
7J F 
^0 F 
50 F 
7J H3L 
10 ^F 
^0 F 
L5 IF 
L^ 112L 
kO IF 
50 111L 
L5 F 
66 IF 
S5 F 

ik ukL 
10 i6f 

50 25L 
26 12^L 
^4-0 F 

LI F 
Yk 111L 
L^ kF 
50 111L 
10 7F 
S5 F 
ho F 
50 111L 
L5 ^F 
10 7F 
S5 F 
it-0 IF 

L5 13L 
Ik h¥ 
h2 k3i 

kl 2F 



NOTES 
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TTT 



u /2 -:i 



N(0) = .35^(1 - U /2) • 2' 



N(l) - (- 1 -fu 1 ) 



1^ 
q_ =, 



N(0) = "ZA 2~ T (n) 



N(l) = H : • 2" T 



Initialize store 
r = N(2) 



LOCATION 


ORDER 


NOTES PAGE 3 


V 1*4- 


35 


50 3F 
79 2F 






36 


Lif 115L 
L*4 F 






37 


Lk F 


(2 n+1 - xr ) . 2"' 






i^-0 2F 


tenrp. store 




38 


50 lllL 

L5 3F 






39 


10 k-F 
66 2F 






ko 


S5 F 








ho ZF 


n-1 




kl 


L5 IF 








LO F 


(N-n)- * 2~7 




k-2 


36 ^-3L — 
22 ^5L 


... . _ 
1 




^3 


L5 2F ^~ 

^0 F 


-J 


store r in (n+l) S3 




kk 


L5 ^3L 

LO ll6L 








^ 


kO k3L 








LQ J? ^ 




k6 


LO 11 5L 
k-0 F 


n = n-1 




*T 


LO 117L : 


n < 1? 






36 51L 


yes 




kQ 


50 lllL 
L5 2F 


no 




k9 


00 IF 








kO 5F 


temp . store 




50 


LL 5? 








36 35L 


if r < 1 
n 




51 


L5 ^3L 
\2 52L 







LOCATION 



52 
53 

5^ 

55 
56 

57 
58 

59 
60 
61 
62 

63 
6^ 

65 
66 

67 
68 

69 



ORDER 



L5 118L 
1+0 F 
L5 F 
kO 5F 
50 111L 
00 2F 
1+0 F 
L5 118L 
kO IF 
L5 2F 
10 k¥ 
kO 2F 
L5 52L 
LO 116L 
1+2 6^-L 
L5 F 
L*J- F 
iA 11% 

14-0 6f 

50 6f 
75 IF 
66 3F 
S5 F 
LO 2F 

ko 6f 

^0 F 
L5 IF 
kO 2F 
L5 6F 
1+0 IF 
L5 6kL 
LO ll6L 
1+0 6hL 
L5 F 
LO 11% 
kO F 



NOTES 
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T~Tk 



3 = 2 



,-7 



N(5) = n . 2 ' at transfer time 



-5 



N(0) = n . 2 



N(l) = j n (2") initially 



N(2) = 3 



n+1 



initially 



Set store address 



(2n + 1) . 2 
temp, store 



-5 



(§*i j ) . 2 - 6 

v x u n' 



temp . store 
(by 59) 



J n+1 J n 



J n J n-1 



decrease store address 



n:n-l 



LOCATION 


ORDER 


NOTES PAGE 5 


V lA 


70 


LO 119L 


n < 






32 5% 


no 




71 


L5 S3 
10 IF 


yes 




72 


kO F 
50 IS 3 


j x 2" 1 — > 




73 


71 3F 
00 k¥ 






7^ 


Lk F 
^0 6f 


(I - x^) 2" 1 




75 


50 120L 
75 3F 






76 


oo 5? 








4o 7F 


$(l) = x/rT mod 2 




77 


U 7F 
50 77L 






78 


26 138L 

40 F 






79 


50 F 






80 


7J 6f 
ko 6f 
L5 7F 


(j_ - x-r ) cos x . 2™ 

j 1 / 




81 


^0 F 

50 81L 


"waste 




82 


26 138L 

kO F 






83 


50 F 








7J S3 


(sin x) j Q . 2" 




8k 


kO F 
50 F 






85 


75 3F 








00 ^F 


(sin x) 3 x , 2 




06 


Lk 6f 








kO IF 


K(l) = ©<. 2" 2 




87 


L5 ^F 
Lk, 121L 


N 





r LOCATION 



ORDER 



NOTES 
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88 

89 
90 
91 
92 

93 
9k 

95 



97 
98 

99 
100 

101 
102 
103 
104 

105 



k2 122L 
k2 123L 
L5 5F 
10 32E 
L^l- 121L 
il-0 5F 
L5 83L 
^2 93L 
k2 95L 
40 F 
50 111L 
L5 F 
10 2F 
66 IF 
S5 F 
^0 F 
F5 93L 
^0 93L 
F5 95L 
ko 95L 
LO 122L 
32 110L 
L5 95L 
LO 5F 
36 93L 
L5 95L 
k2 105L 
k2 10^-L 
LO ll6L 
k2 103L 
^0 F 
50 F 
kO F 
75 F 
00 2F 
^0 F 



end constant 



waste 



J k . 2 



-2 



store J 



end 



waste 



waste 



r n J n-1 °n 



LOCATION 



ORDER 



NOTES 
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106 
107 
108 
109 
110 
111 
112 

113 
21k 

115 
116 

117 
118 

119 
120 
121 
122 



F5 103L 

ko 1Q3L 

F5 lO^L 

kO 10kL 

F5 105L 

kO 105L 

LO 123L 

32 103L 

kO F 

22 F 

00 F 

00 F 

80 F 

00 F 

00 F 

00 3500 0000 OOOOJ 

00 F 

00 62 5000 OOOOJ 

01 F 
00 F 
00 F 
00 IF 
82 F 
00 F 

02 F 
00 F 
Ok F 
00 F 
00 F 

00 3183 0988 6iehj 

35 F 

ko is 3 
S5 F 
kO F 



increment address 



waste 

exit 



-1 



2-7 

2 -39 

(-) 2 

2" 6 

2-5 

i/ir 



,-7 



r 



LOCATION 



123 
±2k 
138 



ORDER 



— ■rUTTW** 



80 2F 
il-0 F 

(S3) 00K 
(T5 )00K 



1 . t^ jf ^m ^ y&A -Stt.il* m 



NOTES 
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V Ik 



.wim,iij"t'»«Am. 



